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Abstract. Extensive calculations in the short-range RVB (Resonating valence bond) subspace on both the 
trimerized and the regular (non-trimerized) Heisenberg model on the kagome lattice show that short-range 
dimer singlets capture the specific low-energy features of both models. In the trimerized case the singlet 
spectrum splits into bands in which the average number of dimers lying on one type of bonds is fixed. These 
results are in good agreement with the mean field solution of an effective model recently introduced. For the 
regular model one gets a continuous, gapless spectrum, in qualitative agreement with exact diagonalization 
results. 

PACS. 75.10.Jm Quantized spin models - 75.40. Cx Static properties (order parameter, static susceptibil- 
ity, heat capacities, critical exponents, etc.) - 75.50.Ee Antiferromagnetics 



1 Introduction 

It is well known that the conventional picture of a long- 
range, ordered, dressed Neel ground state (GS) can col- 
lapse for low dimensional frustrated antiferromagnets. The 
GS of several spin 1/2 strongly frustrated systems has no 
long range antiferromagnetic order and is separated from 
the first magnetic (S — 1) excitations by a gap. The first 
example of such a behavior was given by the zigzag chain 
at the Majumdar-Ghosh point jjj (J2/J1 — 1/2) in which 
case the two-fold degenerate GS is a product of singlets 
built on the strong bonds. 

In some cases the consequences of frustration on the 
structure of the spectrum can be even more dramatic. It is 
now firmly established by many numerical studies that the 
singlet-triplet gap of the Heisenberg model on the kagome 
lattice is filled with an exponential number of singlet states 
||. This property is actually not specific to the kagome 
antiferromagnet (KAF) and could be a generic feature of 
strongly frustrated magnets: It is suspected to occur also 
for the Heisenberg model on the pyrochlore lattice |J , and 
it has been explicitly proved for a one-dimensional system 
of coupled tetrahedra which can be seen as a ID analog 
of pyrochlore Q . 

Since the particular low-temperature dependence of 
many physical quantities is directly connected to the struc- 
ture of this non-magnetic part of the spectrum, many re- 
cent works P,p,|^[^|q,p|, p!o| , pl] | were devoted to understand 
the nature of the disordered GS and low-lying excitations. 
Unfortunately, it is still hard to come up with a clear pic- 
ture of the low-energy sector of the KAF. Resonating Va- 
lence Bonds (RVB) states, for which wave functions are 



products of pair singlets, seem to be a natural framework 
to describe this exponential proliferation of singlet states. 
RVB states were first proposed to describe a disordered 
spin liquid phase by Fazekas and Anderson jl2| for the 
triangular lattice and was reintroduced by Anderson [ fl3| 
in the context of high-T c superconductivity. 

For the kagome lattice, the absence of long range corre- 
lation may lead to consider only Short Range RVB states 
(SRRVB) , i.e. first neighbor coverings of the lattice with 
dimers. The first difficulty which occurs is that the number 
of SRRVB states of a N site kagome lattice with periodic 
boundary conditions is 2 1+(Ar/3) ~ 2(1.26)^ @] whereas 
the number of singlet states before the first triplet of the 
KAF scales like 1.15 N Of course, this does not nec- 
essarily disqualify the SRRVB description but raises the 
question of the selection of the relevant states. 

At the mean field level an answer to this question has 
been given in a recent paper starting from a trimerized 
U) version of the KAF (see Fig. g): 



(1) 



When considering low-energy excitations one can work 
in the subspace where the total spin of each strong bond 
triangle is 1/2. Since there are two ways to build a spin 1/2 
with three spins 1/2, these triangles have two spin 1/2-like 
degrees of freedom : The total spin a, and the chirality r. 
This representation does not simplify the problem because 
spin and chirality are coupled in the Hamiltonian but it is 
no longer the case in the mean field approximation and it 
is possible to solve the mean- field equations exactly 
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Fig. 1. Trimerized kagome lattice: strong bond triangles form 
a (N/3) site triangular lattice. 




Fig. 2. Typical dimer covering of the lattice with nd G f(JA) = 1 
(defaults are marked with a star). 



Low-energy states are SRRVB states on the triangular 
lattice formed by strong bond triangles and their number 
grows like the number of dimer coverings of a N/3 site 
triangular lattice, 1.15 , as can be shown using standard 
methods 

This result was established under the assumptions that 
Jv/Ja is small (trimerized limit) and that quantum fluc- 
tuations can be treated at the lowest order (mean field ap- 
proximation). Therefore two questions remain open: What 
happens beyond mean field approximation ? Can SRRVB 
state give a good description of the energy spectrum in 
the isotropic limit? 

To answer these questions we have studied the KAF 
Hamiltonian in the subspace of SRRVB states with no 
simplifying approximation concerning the non orthogonal- 
ity of this basis. In this subspace the complete spectrum 
is obtained up to 36-site clusters in both trimerized and 
isotropic limit. 

The text is organized as follows: In the first part we 
study the trimerized model and show that mean field pre- 
dictions are robust with respect to quantum fluctuations. 
In the trimerized limit, the low-energy spectrum splits into 
bands in which the average number of dimers lying on one 
type of bonds is fixed and the size of the lowest band scales 
as 1.15^. 

Next we present the results obtained in the isotropic 
limit. Contrary to what was suggested by previous studies 
H the singlet spectrum obtained with SRRVB states is a 
continuum. Moreover the number of states below a given 
total energy increases exponentially for all energy with the 
size of the system considered. 

Finally, we compare the results obtained for KAF with 
the results obtained using the same basis for a non-frustrated 
antiferromagnet, the Heisenberg model on a square lattice, 
and we emphasize the ability of SRRVB states to capture 
the specific low energy physics of frustrated magnets. 

Most of the results presented here contrast with the 
commonly admitted point of view that SRRVB states do 
not provide a good variational basis for this problem. In 



fact, SRRVB states lead to specific numerical difficulties 
due to the fact that they are not orthogonal to each other. 
A way to get around this difficulty is to neglect overlap be- 
tween states under a given threshold. However reasonable 
this approximation may seem, it appears to modify the 
results significantly. It turns out that this approximation 
is not necessary to perform exact numerical simulations, 
even for large systems. In order to clarify this point, some 
technical details about the method we used to implement 
symmetries of the problem and achieve the calculations in 
this non-orthogonal basis are given in an Appendix. 



2 The trimerized model 

As stated in the introduction, the main question with the 
trimerized model (J v / ' Ja < 1) is to know if the mean-field 
selection mechanism (pairing of strong bond triangles) of 
low-lying singlet states is robust when quantum fluctua- 
tions are taken into account. 

Fully trimerized limit - Ground state. Let us start with 
the limit Jv/Ja = 0. In this limit the system consists 
of A/3 independent triangles and the SRRVB GS is ob- 
tained by putting one dimer on each of these triangles. 
Since J v /Ja = this state can be completed to a SR- 
RVB state (A/2 dimers) by putting the iV/6 remaining 
dimer on the J v bonds. The energy of such a state is 
-(3/4)J A (A/3) = -(A/4)J A . In this limit the GS is thus 
obtained by maximizing the number of dimers on the J A 
bonds (A/3). 

By a simple counting argument it is easy to see that 
every SRRVB state contains A/6 = N t /4 triangles, called 
defaults, for which none of the bonds is occupied by a 
dimer (N t — (2 A/3) is the number of triangles): a SR- 
RVB state being a set of N/2 dimer, it leaves (2A/3) — 
(N/2) = (A/6) triangles unoccupied. The number of de- 
faults ndcf(JA) on the J A bonds can take all the values 
from to A/6. In terms of defaults the GS discussed above 
is a SRRVB state which minimizes n^ e {(J^). 
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Fig. 3. The number of coverings with dimers (A) of the 
N/3 site triangular lattice of strong bond triangles is equal 
to the number of covering of the original lattice which satisfy 
ndei{JA) = (sets of (1,2,3) dimer patterns). 
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Fig. 4. In the most general case, when a covering includes 
at least one non-zero bond triangle default, SRRVB states are 
not eigenstates of TL: off-diagonal terms, which overlap with all 
SRRVB coverings, are generated. 



Let us turn to the question of the degeneracy of this 
GS and show that the number of dimer coverings of the 
kagome lattice with ndef(^A) = is exactly the number of 
dimer coverings of the N/3 site triangular lattice formed 
by A triangles. To prove this, we have to check that one 
can associate each GS configuration to a unique dimer 
covering of the triangular super-lattice and vice versa (see 
Fig-J). 

Clearly, to each pairing A of A triangles one can asso- 
ciate a set of dimers (1,2,3) on the kagome lattice. Doing 
so, the number of dimers on A triangles is N/3, which 
is the maximum, and ridef(JA) — 0. Consider now a SR- 
RVB with n^ e f(J A ) = 0. Let us show that there exists a 
unique way to pair A triangles according to the (1,2,3) 
pattern. Starting from dimer 1 on triangle Tl, the exis- 
tence of dimer 2 is necessary because the state is SRRVB 
and the triangle T2 contains dimer 3 because there is no 
default on A triangles by assumption. 

For the triangular lattice, the number of coverings in- 
creases with the number of sites N like Aa^ with at — 
expjj^. f n * f n * ln(4+4 sin x sin y+A sin 2 y)dxdy} ~ 1.5351 
and A ~ 2 |17| . Thus the number of dimer coverings 
of the kagome lattice with Hdef(^A) = increases like 
{al /3 ) N ~ 1.1536*. 

This degeneracy has been obtained considering only 
SRRVB subspace. In the full S = subspace the GS is 
much more degenerate. The model, when J v / J a = 0, sim- 
ply reads: 

W=(J A /2)£{S Ai (S Ai +l)-(9/4)}, (2) 

i 

where S* Ai is the total spin of the triangle i. 

The GS is thus obtained by setting the total spin of 
each A triangle to 1/2 and to couple all the N/3 spin 
1/2 triangles to a total spin of and the degeneracy is 
2^/3 (iV/3)!/[(JV/6)!(l+iV/6)!]. The combinatory factor is 
the size of the singlet sector of N/3 spin 1/2 and the other 
factor refers to the fact that on each of the A triangles 
there are 2 independent ways to build a total spin of 1/2. 



Thus asymptotically the full singlet degeneracy increases 
like 2 2N / 3 /N 3 / 2 ~ 1.5874 N /N 3 / 2 . The table @ summarizes 
the various degeneracies. 



# of sites 


12 


24 


36 


N 


S=0 


132 


2.1 10 5 


4.8 10 s 


2" 
~ F7! 


GS deg. (S=0) 


32 


3584 


5.4 10 5 


1.5874™ 
~ iV3/2 


SRRVB 


32 


512 


8192 


~ 1.26" 


GS deg. (SRRVB) 


12 


72 


348 


~ 1.1536" 



Table 1. Number of singlet and SRRVB states and degeneracy 
of the GS in each of these subspaces as a function of the number 
of sites. 



Fully trimerized limit - Excited states. The situation 
of excited states in the SRRVB subspace is different, even 
when J V /J A = 0, because SRRVB states with n do f(J A ) ^ 
are not eigenvectors of H. In fact this situation occurs 
each time a state includes a default on a triangle with 
non-zero bonds (see Fig. |]). 

Nevertheless, let us consider the results obtained for 
Jv/Ja — (Fig. H). The spectrum splits into bands: the 
first, of zero width, is the degenerate GS discussed above, 
and the other bands consist of linear combinations of SR- 
RVB states with fixed ridei(JA) (a numerical characteriza- 
tion of the dimer coverings in each band is given below). 
The center of each of these bands is — (3/4)7V A J A with 
N A the number of dimers built on A triangles. Since N A = 
{N/3) — ndef(^A), the energy of the center of the 1 + (N/6) 
bands are -(JV/4)J A ,((3/4) - (JV/4))J Al . . . ,-(iV/8)J A . 

Strong trimerization ( J v / J A -C 1). When it is switched 
on, J v acts as a perturbation on the previous spectrum: 
bands with Hdef(>^v) 7^ begins to get wider and to mix. 
In contrast, because it is degenerate when Jv/J A = 0, 
the lowest band is expected to mix with the other bands 
for larger values of J v / J a ■ Let us test this scenario on 
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Fig. 5. Number of states (y-axis) with an energy per site 
smaller than E/N (x-axis) for a 36 site cluster with Jv / Ja = 0. 
The spectrum splits into bands in which nd e f(J&) is fixed. 
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Fig. 6. (a) Same as Fig. q for a 36 site cluster with J v /Ja = 
0.1. (b) Density of state (derivative of (a)). The band structure 
of the spectrum is conserved for a weak trimerization. While 
the excited states band begin to mix the GS band belongs 
separated from the others. 



numerical results for a weak trimerization of the lattice 
(Jv/Ja = 0.1) on a 36 site cluster. Figure || shows the 
spectrum (number of states below a given energy per site) 
and the density of states (DOS) . The DOS exhibit a band 
structure and, as expected, even for such a small value of 
Jv/Ja: gaps between bands are nearly closed. Neverthe- 
less, a very narrow band of states remains very clearly 
separated from the others. 

The existence of this low energy band splited from the 
rest of the SRRVB spectrum indicates that for small values 
of Jv/ Ja i the selection criterion of dimer covering configu- 



Fig. 7. Extrapolation of a(A) for Jv/Ja = 0.1 from 
12,18,24,30,36 site spectra. a(A) exhibit a plateau (1.15) 
which corresponds to the scaling of the lowest band of the 
spectrum. 



rations is the same as for Jv/Ja — 0: the states in the low 
energy part of the spectrum minimize ridef(JA)- I n order 
to test more precisely this scenario let us characterize nu- 
merically the scaling of the bands and verify that /MoK^a) 
is fixed in each band. 

We performed a finite size analysis including all kagome 
clusters with an even number of sites up to 36 sites (12, 
18, 24, 30 , 36). We denote byAf N (A) the number of states 
on a N-site cluster with a total energy smaller than A. For 
all A, the analysis shows that A/at(^) grows exponentially 
with N : 

M N (A) = A(A)a(A) N . (3) 

In the large A limit, since all the states have an en- 
ergy smaller than A, the values of A and a are known 
to be respectively 2 and 2 1 / 3 . Between each band of the 
spectrum, no state comes to increase Af when A increases 
and therefore plateaus appear in a (see Fig. Q). The first 
plateau corresponds to a ~ 1.15, a numerical confirmation 
of what was announced at the beginning of the section. 

Let us turn now to the question of the nature of the 
states in each band. We denote by and N A the opera- 
tors that count for a SRRVB state the number of dimers 
lying on J v and J A bonds respectively. Since on a N site 
cluster, each SRRVB state is made of N/2 dimers, we 
have V v + N A = N/2. Fig. | shows the values of (V v ) 
and (N A ) for each eigenstate of a 36 site cluster from 
the GS to the most excited state. The results are quite 
clear: each band of the spectrum is characterized by a 
fixed value of (N A ) (or (iV v )) which is equivalent to fix 

(n def (JA)) = (N/3) - {N A }. 

SRRVB spectrum versus Exact spectrum. SRRVB states 
on the trimerized kagome lattice spontaneously selects a 
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Fig. 8. Average values of iV A and iV v (y-axis) for each eigen- 
state of TL (from GS to most excited state, x-axis) for a 36 site 
cluster with J v / J a =0.1. 
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Fig. 9. Integrated d.o.s of 12,18,24,30,36 site clusters at the 
isotropic point. The spectrum has the structure of a continuum. 



small set of wave functions (see table [jj) among those 
which minimize energy for Jy/J& = 0. Moreover the num- 
ber of these states scale as the number of singlets in the 
singlet-triplet gap of the KAF at the isotropic limit. If this 
selection is actually relevant, one should be able to identify 
in the exact spectrum at least for a strong trimerization 
the existence of a similar selection. 

To test this point we compare the exact and SRRVB 
spectra for J v /Ja — 0.25 (see table |[ energy per site 
for the 10 first states). The conclusion of this comparison 
is quite clear: The SRRVB subspace reproduces the low- 
energy part of the singlet spectrum and the structure of 
the spectrum (order and degeneracy of levels) is also well 
described. 



Exact Diag. 


SRRVB 


Deg. 


-0.29530 


-0.29513 


1 


-0.29049 


-0.28714 


2 


-0.29027 


-0.28644 


1 


-0.28597 


-0.28453 


3 


-0.28187 


-0.28053 


3 



Table 2. Comparison between exact and SRRVB low energy 
levels for a 12 site cluster with J v / Ja = 0.25. 



In conclusion, beyond mean field approximation, the 
low energy physics of the trimerized KAF is well captured 
by SRRVB states: Low lying states are selected on an en- 
ergy criterion, the maximization of the number of dimers 
on strong bonds, which is equivalent for a weak trimeriza- 
tion to minimize the number of defaults on strong-bond 
triangles. These selected states form a band which con- 
tains a number of states that increases like 1.15^ in agree- 
ment both with exact results and with the effective Hamil- 
tonian approach. 



3 The isotropic model 

When J v / Ja increases to the isotropic limit one may ask 
at least two questions: Does the mechanism described above 
remain valid? Do SRRVB states still provide a good de- 
scription of the singlet sector? 

To answer the first question, we have computed, at 
the isotropic point, the values of (N v ) and (N&) for all 
the eigenstates. The behavior of these quantities is very 
different from the trimerized case: no mechanism tends to 
favor one type of triangle and (iV v ) = (Na) = N/4. This 
means that the simple picture obtained with the trimer- 
ized model is no longer valid in the isotropic case. The 
computation of the spectrum for all even sizes up to 36 
sites confirms the qualitative differences between trimer- 
ized and isotropic model (see Fig. [9]). The mixing of the 
bands which starts for J v < Ja is complete for J v = Ja, 
the band structure has completely disappeared and the 
spectrum is a continuum. 

It is important to emphasize that this result contrasts 
with those obtained for the same model by Zeng and 
Elser[||, who concluded to the presence of a gap inside the 
singlet spectrum. This study was based on an expansion 
using as a small parameter the overlap between SRRVB 
states: the non orthogonality between dimer covering \tpi) 
was neglected under a given threshold of (<fi\<Pj). On the 
contrary, the results presented here involve no approxima- 
tion: the non orthogonality of the basis is fully taken into 
account (see appendix for details). We suspect that the 
difference comes from this approximation. As could be ex- 
pected, our treatment provides a smaller variational value 
of the GS energy. For a 36 site cluster E GS /J = -0.42182 
which is 3% above the exact one (the number of SRRVB 
states is ~ 1.71 1CV 5 of the total singlet subspace). 

The strongest indication that SRRVB states give a cor- 
rect description of the low-lying singlets of the KAF is in- 
deed the continuum structure of the spectrum. Moreover, 
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Fig. 10. a(A) (see Eq.jj) for the standard kagome model. In- 
set: zoom of the circled region. The abscissa of the points de- 
noted by the crosses are the exact diagonalization values of the 
singlet-triplet gap for the 36 site and 12 site clusters. 



the shape of the spectrum is very similar to the one ob- 
tained by exact diagonalization (ED) II. In order to test 
this point more precisely, we have again computed a(A) 
at the isotropic point 



(see Fig. [10 ) . Plateaus no longer 
appear in a(A), which confirms the complete mixing of 
the bands. More interestingly, this analysis shows that, 
for all A, the number of SRRVB excitations increases ex- 
ponentially with the size of the systems. This proves that 
SRRVB states not only reproduce the continuum nature 
of the spectrum but give a good description of the expo- 
nential proliferation of singlets states in the low energy 
sector of the KAF. 

Since the SRRVB subspace cannot give information 
about magnetic excitations the question of the counting 
of states below the first triplet is rather delicate. To dis- 
cuss this point, one has to take the exact singlet-triplet 
gap value to make the counting in the variational SRRVB 
spectrum. Doing so on has to keep in mind that even if the 
SRRVB spectrum gives a good description of the structure 
of the low lying singlets (order of levels, degeneracy, ex- 
ponential proliferation) the energy scale of the excitations 
above the GS might be slightly different from the exact 
one: SRRVB arc not the exact eigenstates of the Hamil- 
tonian which are more probably dressed SRRVB states 
including fluctuations that modify the energy scale. But, 
given the relative accuracy of the GS this point should 
not prevent us from doing a semi-quantitative comparison 
between exact and SRRVB results. 

In fact, for a 12 site cluster we have checked that the 
low energy structure of excitation spectrum is correct for 
J V /J A = 1 up to the first triplet state (see Table ||) . 

For a more general quantitative discussion on the pro- 
liferation of low-energy singlets let us analyze the shape of 
a(A) (see Fig. [h]) obtained from SRRVB spectra of 12, 18, 
24, 30 and 36 site clusters. The range of the exact singlet- 
triplet gap extends from ~ 0.38 J for the 12 site cluster to 
~ 0.17 J for 36 sites {§ which corresponds to the circled 
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Table 3. Comparison between exact and SRRVB low energy 
levels for a 12 site cluster at the isotropic point Jv/Ja = !• 
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Fig. 11. Integrated d.o.s. for the Heisenberg antiferromagnet 
on the square lattice in the SRRVB subspace for 16, 18, 20 and 
26 site clusters. The low energy structure of the spectrum is 
no longer a continuum. 



region and the inset of Fig. |lC]. It is remarkable that in 
this energy range the value of a for SRRVB spectra goes 
from ~ 1.18 to ~ 1.15 which in good agreement with ED 
scalings. 



4 Discussion 

At this point, it is fair to ask whether the continuum struc- 
ture of the spectrum obtained with SRRVB states is re- 
ally a specific feature of frustration captured by this basis 
or simply a generic characteristic of the spectra that such 
states would provide on any lattice. To answer this impor- 
tant question let us compare the results on the kagome lat- 
tice with the SRRVB spectrum for a non-frustrated model, 
the Heisenberg model on the square lattice (see Fig. |i~l| ). 

The structure of the SRRVB excitations on the square 
lattice is qualitatively different from the structure obtained 
for the kagome lattice: In particular there is a gap between 
the singlet GS and the first excitation. Even if it is seems 
difficult to extract a precise value for this gap (see Fig. [l2|) , 
the finite size analysis strongly suggests that it remains fi- 
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Fig. 12. Total energy difference between SRRVB first excited 
state and GS for the square lattice (16, 18, 20 and 26 site 
clusters) . The Data suggest a non-zero value for the gap in the 
thermodynamic limit. 



nite in the thermodynamic limit. Of course, this does not 
describe the actual singlet spectrum of the square lattice, 
which is gapless due to two-magnon excitations. But it 
shows that the structure of the RVB spectrum is specific 
to very frustrated lattices. 

In conclusion, SRRVB states on the kagome lattice al- 
low to capture the specific low energy properties of the 
model in both trimerized and isotropic limits. In the trimer- 
ized model it gives a simple picture of the non magnetic ex- 
citations and a selection criterion of the low-energy states 
which are built by minimizing the number of defaults on 
strong bond triangles. The number of such states increases 
like 1.15^. The states matching this criterion can also 
be seen as short-range dimer coverings of the triangular 
lattice formed by strong-bond triangles which confirms, 
beyond mean field approximation, the relevance of the ef- 
fective model approach. At the isotropic point, SRRVB 
states lead to a continuum of non-magnetic excitations in 
accordance with ED results. Moreover the shape of the 
SRRVB spectrum is very similar to the exact one and the 
number of low-lying singlets increases exponentially for all 
energy range with the size of the system considered. 

Finally, these properties of the SRRVB spectrum are 
not just a property of this kind of states since the SRRVB 
spectrum has a gap in the case of the square lattice. So 
one may conjecture that they provide a good description 
of the low-energy singlet sector of very frustrated magnets 
only. Work is in progress to test this idea on the pyrochlore 
lattice. 
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5 Appendix: numerical method 

Working with SRRVB states as a truncated basis leads to 
non-trivial numerical difficulties which, as paradoxical as 
it may seem, make the problem of the determination of the 
spectrum more tricky in this truncated subspace than in 
the full space of spin configurations. This is a consequence 
of the non-orthogonality of RVB states. 

If and \(fij) are two SRRVB states, the overlap is 
given by fill 



(4) 



where rib((pi, tfj) is the number of closed loops in the dia- 
gram where the two states are superimposed and s(</?j ,tpj) = 
(— l) p , where p is the number of misorientcd dimers com- 
pared with the reference orientation (see Fig. |l3| ). 

In the case of a non-orthogonal basis, the eigenvalues 
are solutions of the so-called generalized eigenvalue prob- 
lem, 



dBt((<p i \H\(p j )-E(<pi\<p j )) = 0, 



(5) 



in which the overlap between states appears explicitly. 
Since we are interested in the structure of the spectrum, 
we need to diagonalizc completely the Hamiltonian and 
therefore iterative techniques (typically Lanczos) must be 
avoided. On the other hand, solving (|E|) with standard 
routines, one is limited to small systems. 

To achieve a complete diagonalization for large systems 
(typically 36 sites) it is crucial to take into account all the 
symmetries of the system in order to break the Hilbert 
space into smaller pieces. This technique is indeed very 
standard but is usually used in a context where the basis is 
orthogonal (e.g. spin configurations) which makes it quite 
convenient. The non-orthogonal case is far less simple and 
is worth paying some attention. 

The aim of this appendix is to explain how it is pos- 
sible, starting from a set of configurations that can be 
non-orthogonal, to build an orthonormal basis of vectors 
that are eigenstates of all the symmetries of the problem 
in each symmetry sector. Since this linear algebra problem 
is planned to be solved numerically one is interested in re- 
ducing as much as possible the information to be handled. 
Therefore one does not work explicitly with this orthonor- 
mal basis but with linear combinations of suitably chosen 
configurations called representatives. 

The text is organized as follow: we define the repre- 
sentatives, we show how the number of representatives 
has to be reduced depending on the symmetry sector and 
finally explain how one can go from representatives to the 
orthogonal basis of the symmetries eigenvectors. 

Representatives. Let us denote by N$ the order of the 
symmetry group of the system and Si, i = 1, . . . , N$ the 
elements of this group. It is possible to make a partition of 
the set containing all the configurations in subsets where 
configurations are related to each other by a symmetry 
Si. Each of these N p subset can be represented by a con- 
figuration \pi), i — 1, . . . , N p , called representative, of the 
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Fig. 13. The overlap between two dimer coverings depends 
on the number of loops in the transition graph (bottom) . Note 
that the relative orientation of dimer coverings in this figure 
corresponds to the reference orientation which produces a pos- 
itive overlap. 



subset since by construction all the others can be obtained 
by applying symmetries on it. From a numerical point of 
view the set of the representatives is the minimal infor- 
mation needed. 

Reduction of the number of representatives in a given 
symmetry sector. In this section we will consider a given 
symmetry sector s characterized by a set of characters 
Xs(Si), ■ ■ ■ ,Xs(Sn s )- We are going to show that it is not 
necessary to keep all the representatives to generate a ba- 
sis in the sector s. 

Let us consider a given linear configuration of repre- 
sentatives, 

N p 

\i,)=Y j a i \ Pi ). (6) 

i=l 

The true state of s associated to \ip) is given by |20|] , 

, N s n p N s 

v S i=l i=l p=l 

where |s p (pj)) stands for the image of the configuration 
Pi by the symmetry S p . For a given representative \pt), 
let us denote by £i the set of indices q of the symmetries 
that leave the configuration pi invariant (s q (pi) = pi), and 
£i the remaining indices. With this notation take the 
form, 

i=l \q££i J 

N p 

+ ^^a, Xs (S,)K(p,)). (8) 

Let us denote by Q s the list of indices i of the rep- 
resentatives pi such as J2q£Ei Xs{Sq) — in the sym- 
metry sector s. It is obvious to note that all the repre- 
sentatives with an index in Q s disappear from the first 
term of eq. (||). What we are going to show is that they 
also disappear from the second one. Let be Pi such as 
Yjq&SiXs{S q ) = and p E £ l . One has Y, qe e z S i\P^) = 
J2g & £ l X*s(S p )S p S q \p l ) = J2 q e£ t X*s(Sp)S q \p l ). This result 
does not depend on p S £, and thus one does not modify 
the result by applying [1/Card(£j)] X) P e£i on * ne P rev i- 
ous expression which proves what was announced. This 
leads to a reduction of the number of representatives in 
the symmetry sector s, namely N s = N p — Card(Q s ). 

From non- orthogonal representatives to the orthonor- 
mal basis of symmetries eigenvectors. In the general case 
of non-orthogonal representatives, it is convenient to intro- 
duce a mixing matrix fi s in order to build the orthonormal 
basis of symmetries eigenvectors. We will considerer from 
now linear combinations of mixed representatives, 



3=1 



(9) 



M. Mambrini, F. Mila: RVB description of the low-energy singlets of the spin 1/2 kagome antiferromagnet 



9 



All the problem is to chose fx s such as the symmetrized 
states \l) of \l) according to (0) form an orthonormal basis: 
(l\rh) — Sim- Let us show how this condition writes in the 
cases of orthogonal and non-orthogonal representatives. 

Orthogonal case. In this simple case where (jpi\pj) = 
Sij, the condition (l\fh) = 5i m is, 

N s Ns 

Sl m = Y Y Vjl " Vmi(Pj\ S q\Pi) 

«,i=i <3 =1 

= Y Pit Id 5 *i 

N s 

= Y M^Mmidegfe)^-, (10) 

where Xs(^p) is the character of the symmetry S p in 
the sector s, deg(pi) the degeneracy of representative pi 
(i.e . the number of symmetries under which it is invariant 
PjJ), s p(Pi) the image of the configuration pi by S p , and 
N s the size of the sector s. 

It is easy to see that 

lAj = Sij/y/degfa) (11) 

fulfills © 

General case. In the non-orthogonal case, the condition 
( ^0| ) now reads, 

Sim = Y I'Y'l'n.'l-* ( 12 ) 

where, 

p=i 

Here again, the indices i and j runs from 1 to N s (the 
size of the sector s). To determine /i s we diagonalize / : 

pt/P = Diag(d 1) ...,d w J (14) 
One can check that, 

/4 = JL Pa (15) 

satisfies condition (p^). 

The basis is orthogonal and in this new basis the 
Hamiltonian is block diagonal, each block corresponding 
to one symmetry sector. Thus, it only remains to diago- 
nalize the Hamiltonian in each of these representations to 
get the whole spectrum : 

N s N s 

(l\H\rh) = Y Y^^XsiSp^Pjlnls^)) (16) 



In conclusion the procedure described above turns the 
generalized eigenvalue problem of n x n matrices into 2 x 
N s conventional diagonalizations of ~ (n/N s ) x (n/N s ) 
matrices. The point we now want to stress is that the 
treatment described above is exact and does not introduce 
approximation. The subspacc of RVB states is a truncated 
subspace in the sense that it is not stable with respect 
to an application of the Hamiltonian and only the RVB 
restricted Hamiltonian is studied. But, the use of symme- 
tries does not act as a new restriction of the Hamiltonian 
in each representation of the symmetry group. In the new 
basis the Hamiltonian and the overlap matrix are actually 
block diagonal, each block corresponding to one represen- 
tation. Thus the spectrum obtained as well as mean values 
of operators calculated are the same as those one would 
obtain by solving with brute-force the generalized eigen- 
value problem if it was possible. 
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